MR spectroscopy with parabolic magnetic field: spin-oscillator coupling effect 
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The spectrum of the spin particle in oscillatory potential 
subjected to external parabolic magnetic field B = (_Bo + 
Gx + Gx^)i is obtained. The structure of energy levels of 
the considered system allows to identify the frequency of the 
oscillator via the spectrum of spin sublevels coming only from 
one oscillatory level. The effect is due to the gradient terms 
in the form of the field. 



In a recent decade one can observe a rapid grow of 
interest in the Stern-Gerlach (SG) interaction ^ within 
MR spectroscopy . The main reasons are the follow- 
ing: (i) partial saturation of the possibilities of the con- 
ventional MR methods based on Bloch-Purcell paradigm 
1^ (low sensitivity, especially for the nuclear MR (NMR) , 
growing costs of aparature, existence of cheaper, com- 
petitive methods like e.g. ultrasonic ones), (ii) recent re- 
ports on new and promising effects (e.g. atomic force mi- 
croscopy 1^,^ , MR microimaging , MR spectroscopy 
for selected Zeeman states (iii) the need for new so- 
lutions in the domain of NMR quantum computing (the 
main problem in the recent implementations of quantum 
computing is low signal intensity) |^ (iv) the recent de- 
velopment in the domain of ferromagnetics and super- 
conductors allowing to produce strong magnetic fields of 
required shapes. In this situation one is seeking for new, 
both theoretical and experimental solutions for MR, in 
particular by exploiting the SG interaction. 

Recently Sidles has proposed the interesting model 
of spin-i particle mounted to the harmonic oscillator. 
The particle is assumed to interact with the magnetic 
field via Stern-Gerlach Hamiltonian plus some time- 
dependent, constant in space magnetic field. He solved 
the model and showed that the wavepacket of the particle 
exhibits the splitting into two parts moving in opposite 
directions. He also predicted the importance of possible 
oscillator based NMR spectroscopy for molecular imaging 
devices. 

In this Letter we consider the coupled spin-oscillator 
system in external heterogeneous magnetic field. The 
model is solved exactly and provides physical spectrum 
in which the oscillator frequency is coupled to the field 
derivatives parameters. It implies that in principle one 
can recognise the frequency of oscillator basing only on 
signals coming form transitions between sublevels of one 
oscillatory level. 



Let us consider the particle that subjects to the sim- 
ple dynamic due to the harmonic oscillator Hamiltonian 
T-L^g^{x — a) . Here a is the position of the minimum of the 
oscillator potential with respect to the origin of the refer- 
ence frame (we imagine that the particle is in a sample of 
finite dimension / situated in the origin of the reference 
frame, so that a can vary within the region (— /, /)). The 
Hamiltonian of the particle 
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has the following spectrum and corresponding eigenvec- 
tors: 
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with Hermite polynomials 

Suppose now that the particle has the spin S and 
recall that the quantum operator corresponding to the 
projection of the spin onto the z axis is described as 
i = X]m=-s ^'^^^ the corresponding spec- 

trum and eigenvectors TiM, \M). 

Consider the inhomogeneous magnetic field 



B{x) =B{x)z= {Ba + Gx- 



Gx^)z 



(3) 



along z axis. Here G, G stands for the gradient and 
second derivative parameters. Note that the field has the 
value Bq = Bqz for a; = an that in general it depends 
quadratically on the spatial coordinate x. 

Now if we put our particle in the field (||) then its 
Hamiltonian becomes 



n = K,,{x 
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In the above formula the difference between spatial and 
spin degrees of freedom has been stressed. In particular 
the identity operator / has been used. The new energy 
spectrum of our particle can be found seeking the eigen- 
vectors in the form ]((>) ® |M) (M = -S, -S + 1, S). 
This leads effectively to the sequence of 2S + 1 M- 
dependent shifted quantum oscillators. Eigenvalues each 
of those oscillators can be easily solved and the final en- 
ergy spectrum of the particle is given by 
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EM,n = nJ - ^^'^^^\ n + i) - 7(^0 + Ga + Ga^)hM 
\ m 2 
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and the corresponding eigenvectors can be written in the 
form 



$M,« = \4'M.n) ® \M) 



(6) 



The spatial coordinate function is defined by the oscilla- 
tor eigenvectors (||) in the following way 

J. r ^ T /o j{G + 2Ga)n^M. 
q)M.n{x) = ^n{^lM;x - a ^- ) (7) 



where ^ " 

It is important to note that the squared component of 
the field represented by its second derivative G must not 
be too large as then it suppresses the squared term of 
the particle oscillator potential ™ ^ leading to Hamil- 
tonian unbounded form below. This occurs for the values 
violating the inequality 



G < 
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where M^ax — S. In such case the first term of the spec- 
trum (||) becomes imaginary for some M. Now, if the 
particle were a light, spin endowed part of molecule (the 
model which we shall discuss subsequently) , this could be 
interpreted as the dissociation of the molecule caused by 
the strong gradient of the magnetic field (herewith, we 
will keep, a bit roughly, the term "dissociation"). The 
last inequality establishes the boundaries of our model. 
In this context it is convenient to introduce the new dis- 
crete parameter defined by 



j^^2jGnM 
Using M the inequality (||) writes as 



< 1 



(9) 
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i.e. the dissociation occurs if M reaches the unity. In 
terms of M, the dissociation occurs if M reaches unity, 
so that the scope of the presented model is described by 
inequality M < 1. Roughly speaking M is equal to the 
ratio of M to the number of the maximal spin level for 
which dissociation occurs. As we are interested in the 
effects caused by gradients, we put Bq = 0. Then the 
energy spectrum can be written in the following form 



EM,n = hQ{- 



.1 



1-M- 



mn^a^ M mn\G/Gy 



mn^{G/G)a M 
2 1 -M 



1 - M 



4(1 - M) 



(11) 



Consider first the case of infinitely small sample con- 
centrated in the origin of the reference frame (i.e. put 
a — 0). Then we obtain especially appealing form of 
energy spectrum 



EM.nia = 0) = E^,^,{n,nWl^M-Ef^^in, G/G) ^^^__ , 



where E^^^ = hn^{ 



2/1 
2 



(12) 

and E^,i,{n,G/G) = 



Here we have two competitive terms: the first one rep- 
resents the energy of quantum oscillator of frequency f2 
(i.e. the energy of the particle in absence of the field) 
while the second one is nothing but the total energy of 
the classical oscillator of the same frequency, of the am- 
plitude determined by the shape of the field. Both terms 
have weights depending on the scaled spin number. Let 
us now keep the value of G constant (suitably chosen, 
in order not to be too close to the dissociation regime, 
e.g. to have M ~ 0.1). Now varying the linear gradient 
G we obtain smooth transition from quantum oscillator 
regime to the classical one. Provided G is sufficiently 
large, so that the classical part dominates, we have a 
kind of amplification of the vibrations of the particle by 
the magnetic field. In both cases the frequency H. de- 
termines the ratio of splitting of the spin levels. Thus 
by measuring the latter, we can obtain the information 
about the value of the frequency. Note that, remarkably, 
the same is impossible in homogeneous field: the possi- 
bility of monitoring vibrational frequency via magnetic 
resonance is exclusively due to the gradient terms. 

Note that the double "quantum-classical oscillator" 
structure of the energy levels is due to quadratic gra- 
dient term. Indeed, for G = we we cannot use the 
scaled spin level number M as the latter is defined for 
G ^ 0. In such a case we have no "dissociation" as there 
is no denominator singularity in formula (|^). If we put, 
instead, G = (with a 7^ we still have two oscillators, 
however the latter one follows from finite dimensions of 
the sample (so that the centre of the oscillations may feel 
nonzero value of magnetic field) 



EmAG = 0) = E^,^,{n, nWl-M- E^J,,{n, a) 
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On Fig. 1 we present these levels for electronic MR 
(EMR) versus quadratic gradient G for the system of spin 
S — 3/2 (we keep G being nonzero, but weak). One can 
observe remarkable crossing points for sufficiently large 
G. Again, for a given value of G the level difference is 
determined by the vibrational frequency fl. 

Let us summarise our results. The very important fea- 
ture of the obtained spectrum of the particle is that the 
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spin quantum number M has been coupled to the charac- 
teristic frequency of the particle. This means that now 
for fixed quantum oscillator number n the M-dependent 
energy sublevels form the structure which depends on the 
characteristic frequency oj the particle. Indeed, putting 
G = G = one can immediately see that it does not 
happen for uniform field B(a;) — Bq. It is remarkable 
result for the following reasons. Imagine the molecules 
with spin-less heavy "core" and spin 5* light part cou- 
pled to it. Assume that we have a mixture of unknown 
molecules of that kind, each of them possessing unknown 
characteristic frequencies Oi, il2, ... . Then the fact that 
the spectrum (|]) (with nonzero G or G) depends on the 
frequency will provide the possibility to identify them 
basing directly on the NMR techniques that generically 
deal with 2S' -I- 1 spin sublevels of the given level. 

It is interesting that if we have G — Q and nonzero 
G then the n dependent term does not depend on quan- 
tum number M . Of course, the structure of spin sub- 
levels does not depend on n. Then all the molecules of 
one kind (distinguished by fi) would give the same signal 
from M sublevel of all oscillatory levels. A very interest- 
ing feature of the considered interaction is the emerging 
amplification of the vibrational mode to the form of clas- 
sical oscillator. The effect is due to nonzero G. The 
amplitude is determined by the shape of the field, and 
the amplification is obtained for high values of the linear 
gradient. 

Concerning the experimental realiation of the pre- 
sented spin-oscillator coupling effect note that for the 
typical molecular frequences (« IQ^^Hz) one would need 
enormously strong gradients to observe the effect. In- 
deed, the strongest gradients obtained so far amount to 
10®T/to Now, for molecule with incorporated atom 
of oxygene the above gradients allow to realize the ef- 
fect for frequences up to \Q'^Hz. Besides the highest 
frequency for which the effect can be observed satisfies 
roughly G ^ Vl^ . On the other hand, there exist molec- 
ular systems that exhibit low frequences Thus the 
effect is difficult but still possible to realize experimen- 
tally. 

In conclusion, we have demonstrated that MR spec- 
troscopy with parabolic magnetic field provides subtle 
spin-oscillator coupling effects, that allow to select the 
characteristic frequency of the system. It seems to open 
new prospects for the MR applications such as quantum 
level separations or quantum computation Fi- 
nally, we expect that the idea of incorporating higher 
order gradients will prove fruitful in looking for new so- 
lutions within MR domain. We also believe that the ob- 
tained results will stimulate development of techniques 
of producing strong gradient field of different shapes. 
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FIG. 1. Structure of the EMR energy levels (a) in parabolic 
field for finite dimension of the sample (5 = 3/2, a = 10~*m, 
n = IQ^Hz, G = -0.003T/m) and oscillatory spectrum (b). 
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